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Abstract We simulate evolution of cometary H II regions based on several champagne flow 
models and bow shock models, and calculate the profiles of the [Ne II] fine-structure line 
at 12.81pm, the H30a recombination line and the [Ne III] fine-structure line at 15.55pm 
for these models at different inclinations of 0°, 30° and 60°. We find that the profiles in the 
bow shock models are generally different from those in the champagne flow models, but the 
profiles in the bow shock with lower stellar velocity (< 5km s“^) are similar to those in the 
champagne flow models. In champagne flow models, both the velocity of peak flux and the 
flux weighted central velocities of all three lines are pointing outward from molecular clouds. 

In bow shock models, the directions of these velocities rely on the speed of stars. They have 
the similar motion in high stellar speed case but opposite directions in low stellar speed case. 

We notice that the line profiles from the slit along the symmetrical axis of the projected 2D 
image of these models are useful for distinguishing bow shock models and champagne flow 
models. It is also confirmed by the calculation that the flux weighted central velocity and the 
line luminosity of the [Ne III] line can be estimated from the [Ne II] line and the H30a line. 
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1 INTRODUCTION 

Massive stars always form in dense molecular clouds and transfer a great amount of energy from their 
ionizing fluxes and stellar winds to the surrounding interstellar medium (ISM). This affects the kinetic 
and thermal energy of the molecular clouds so as to change the velocity structure and the morphologies 
of the clouds substantially. The ionizing photons from a massive star can ionize the surrounding ISM and 
form an H II region. In a uniform environment, the H II region will be spherical and is called Stromgren 
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sphere (Stromgren 1939). However, observations show that a large number of young H II regions have 
a cometary morphology (Wood & Churchwell 1989; Kurtz et al. 1994; Walsh et al. 1998). Tenorio-Tagle 
(1979) pointed out that these cometary HII regions result from the density gradients in the molecular cloud. 
In this explanation, the comet-shaped H II regions are called blister H II regions or named as champagne 
flows. The Orion Nebula is the archetype of this kind of H II regions. Observation shows that there is 
a bright ionization front on the surface of a molecular cloud and ionized gas flows away from the cloud 
(Israel 1978). The champagne flow model was first created by Tenorio-Tagle and coworkers (Tenorio-Tagle 
1979; Bodenheimer et al. 1979; Yorke et al. 1983). In these studies, the H II region is assumed to form 
in a uniform molecular cloud, but close to the boundary of it. The H II region will expand and break out 
into the intercloud medium of low density to form a shock. The Mach numbers of the resulting shocks are 
shown to be 30 — 42, and the velocities of the ionized gas can reach up to 40 km Since then, a series 
models were developed by including the effect of a stellar wind, a magnetic held or an exponential density 
distribution (Comeron 1997; Arthur & Hoare 2006; Gendelev & Krumholz 2012). In all of these models, a 
density gradient is assumed to exist, which is considered as the main cause of the cometary morphology. 

There is an alternative explanation to the cometary shape of H II regions called bow shock. In this kind 
of models, a wind-blowing ionizing star moves supersonically with respect to the dense molecular cloud, 
and forms a shock in front of the star. A cometary H II region can also form in this case (Mac Low et 
al. 1991; van Buren & Mac Low 1992). Wilkin (1996) created an analytic model to derive the bow shock 
characteristics. Other researchers made a series numerical models to simulate bow shocks (Comeron & 
Kaper 1998; Arthur & Hoare 2006). The presence of the density gradient seems unavoidable in the non- 
homogeneous environment in molecular clouds. High speed stars are not common in the Galaxy. But, it can 
happen in the region with a high stellar density, such as massive star forming regions. Hence, it is necessary 
to compare the two possible causes of cometary H II region to interpret observations. 

Arthur & Hoare (2006) pointed out that simple champagne flows without a stellar wind will not show 
a limb-brightened morphology. By contrast, bow shock models always display this morphology. They also 
pointed out that the two kinds of models can be distinguished in the kinematics. For example, the highest 
velocity of ionized gas with respect to the molecular clouds is at the head of the cometary H II region 
in the bow shock models while it exists at the tail in a champagne flow model. It is the most notable 
that the directions of motion of the ionized gas at the head are always toward the molecular clouds in 
bow shock models but could be away from the clouds in champagne flow models (Arthur & Hoare 2006). 
Since the profiles of emission lines are partly determined by kinematics, this difference in kinematics could 
lead to different line profiles. Line profiles are direct observing properties that can be accessed with high 
resolution observations. Therefore, line profiles can be used as the criterion to select better models for 
certain observations. 

In this paper, we present the profiles of the [Ne II] 12.81pm line, the [Ne III] 15.55pm line and the 
hydrogen H30a (31 — 30) recombination line from cometary H II regions by simulating bow shock models 
and champagne flow models with different parameters. The [Ne II] 12.81pm line is the brightest line from 
HII regions in the mid-infrared, which can be accessed by TEXES on IRTE or EXES on Stratospheric 
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ground. In order to investigate the profile of the [Ne III] 15.55pm line, we have to rely on future space 
telescopes. The H30a line can be observed by Caltech Submillimeter Observatory (CSO). These lines 
are useful to study the gas kinematics of compact HII regions. We compare the line profiles between the 
bow shock models and champagne flow models and discuss the possibility to distinguish these two kinds of 
models based on these line profiles. The organization of this paper is as follows; In Section2, we describe the 
method of our numerical simulation. In Section3, we present the results of the numerical models. Section4 
show our conclusions. 

2 METHOD 

2.1 Method of Hydrodynamics, Radiative Transfer and Thermal Processes 

It is necessary to simulate the time evolution of the cometary H II regions in order to obtain the profiles of 
lines of interest. To correctly describe the time evolution of the regions, one has to treat the transportation of 
energy properly. In the situation of emission nebula, the gravity is not important, but the hydrodynamics and 
the radiative transfer are essential processes for energy transportation. In this paper, both hydrodynamics 
and radiative transfer are considered. 

A 2D explicit Eulerian hydrodynamic method is used to treat the evolution on a cylindrically symmetric 
grid. Most of the models in this work are computed on a 250 x 500 grid. A big grid of 400 x 400 is used 
when computing blister H II region model to check the contribution of the ionized gas far from the star at 
the sides to the emission line profile. The results show that the grid of 250 x 500 is enough. The cell size 
of the grids is chosen to be dr = 0.005 pc. A HLLC Riemann solver (Miyoshi & Kusano 2005) is used to 
solve the hydrodynamic conservation equations. 

When treating the radiative transfer, we consider the star as a single source for ionizing and dissociating 
radiation. We solve the radiative transfer for EUV (hv > 13.6eE) and EUV (11.26ey < hv < IS.GeE), 
respectively (Diaz-Miller et al. 1998). The black body spectrum is assumed for the ionizing star. We as¬ 
sumed the ”on the spot” approximation when treating ionizing photons. Eor dissociating radiation, because 
the column density of molecular hydrogen exceeds in our models, the EUV lines are optically 

thick. Thus, the self-shielding by H 2 becomes important (Hollenbach & Tielens 1999). The dissociation 
rate and reformation rate of hydrogen molecules are calculated by using a simple self-shielding approxima¬ 
tion introduced in Draine & Bertoldi (1996) and the method in Hollenbach & Tielens (1999), respectively. 
The dissociation and reformation of CO molecules are also included in the models, following the methods 
given by Lee et al. (1996) and Nelson & Langer (1997). 

In the ionized region, the photoionization heating is considered as the only heating process (Spitzer 
1978). We use the cooling curve for solar abundances given by Mellema & Lundqvist (2002) to compute 
the radiative cooling rate. This cooling curve is derived on the assumption that the cooling of the gas is 
due to collisional excitation of hydrogen and metal lines and hydrogen recombination. Outside of the H 
II region, the heating processes of the gas include photoelectric heating, heating from photodissociation, 
reformation of hydrogen, cosmic ray and EUV pumping of H 2 molecules as heating processes. The fol¬ 
lowing cooling processes are considered; atomic fine-structure lines of [O I] 63p,m, [O I] 146p,m and [C 
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lisions (Hollenbach & McKee 1979; Tielens 1985; Hollenbach & McKee 1989; Bakes & Tielens 1994; 
Hosokawa & Inutsuka 2006). Since the results of the photodissociation region will not be presented in this 
paper, the purpose of including the radiative transfer of the photodissociation radiation is mainly to conserve 
the energy and momentum of the gas. 

2.2 Line profiles 

After we solve the continuum, momentum and energy equations of the hydrodynamic models, the method 
derived by Glassgold et al. (2007) is applied to compute the line luminosity for the [Ne II], [Ne III] and 
H30a lines. The line luminosity i at a given velocity v is as follow; 


= < 


^k^k^ulhl^ul 

AbNeX{Ne+)nPuAuihiyui 

AbNeX{Ne'^+)nPuAuihvui 


for H30a, fc = 31, * = 0 
for [Ne II], i=l 
for [Ne III], i = 2 


( 2 ) 


where vios{r) is the line-of-sight component of the velocity vector, hi/ui is the photon energy, vth is the 
thermal velocity of gas. b^Nk is the number density of hydrogen atoms in the kth electronic energy level 
{k = 31 for H30a). Nk is the theoretical value for the number density in the kth level expected in LTE and 
is proportional to nf,. bk is departure coefficient (Seaton 1959). n is the number density of H nuclei. The 
abundance of Ne, AbNe = 1.0 x 10“^, is adopted in these models (Holweger 2001). Aui is the Einstein 
emission coefficient for corresponding transitions (Alexander 2008; Glassgold et al. 2007). X{Ne~^) and 
X{Ne‘^~^) are the fraction of and ions, respectively. Pu is the excitation fraction of the upper 

state and is computed as in Glassgold et al. (2007); 


j [2(1 + n^r/ne)exp{ll22.%/T) + 1 ]-^ for [Ne II] 

Pu — \ 

[ [1 + (5/3)(l + ncrilne)exp{925.3lT) + (1/3)(1 + ncr2/ne)exp{-399/T)]-^ for [Ne III] 

(3) 

where the critical density Ucr = 5.53 x 10^T°-® cm~^, Ucri = 3.94 x 10^T°'^ cm~^ and ncr 2 = 7.2 x 
^Q2y0.5 relative fractions of Ne, Ne'^ and iX{Ne) + X{Ne~^) + X{Ne'^~^) = 1) 

are computed through the ionization-recombination balance equations. The photoionzation cross sections, 
the recombination rate coefficients and the charge exchange rate coefficient are given by Henry (1970), 
Pequignot et al. (1991) and Glassgold et al. (2007), respectively. The collisional ionization of Ne is not 
considered in photoionized region where temperature is much lower than the critical temperature Tc = 
2.50 X 10® K, and the gas density is relatively low. The contribution of the hot stellar wind bubble (T = 
10®“® K) to line luminosities is negligible because of its low density (n < 5 cm“®). Hence, collisional 
ionization of neon atoms is not considered. In addition, since the average energy of photons from young stars 
are not high energy, the fractions of highly ionized ions as Ne^~^ and A^e"‘“'' are generally not considered 
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the relative fractions of ions and other higher ionization state Ne ions. They are always lower than 

0.002 in photoionized region, and can be safely neglected. X rays and EUV photons with energies greater 
than 21.56eV can both photoionize neon, but we neglect X-ray ionization because of the lack of x ray 
sources relative to EUV photons. The results of the H30a recombination line can be applied to other H 
recombination lines. They have the same normalized profile when the pressure broadening effect can be 
neglected. 

In this paper, we assume that the observers are viewing the cometary regions from the tail to the head. 
Therefore, a blue-shifted velocity suggests that the gas mainly moves at a direction from the head to the tail 
and a red-shifted velocity suggests the gas mainly moves at the opposite direction. 

3 RESULTS 

In this section, seven models are presented. Eour of them are bow shock models, and the rest are champagne 
flow models. The parameters of these models are presented in Table 1. These models are selected to test 
the effects of different model ingredients on the line profiles. The values of stellar parameters are based 
on Diaz-Miller et al. (1998) and Dale et al. (2013). The initial density is consistent with the condition in 
compact H II region. In the following section, the models will be described and analyzed individually. 


Model 

M{MQyr-^) 

Vi,{km s ^) 

Vw {km s 

log[Sc/v(s-l)] 

no (cm 

Scale Height(pc) 

A 

9.93 X 10"'^ 

10 

2720.1 

48.78 

8000 

0 

B 

9.93 X 10"'^ 

0 

2720.1 

48.78 

8000 

0.05 

C 

9.93 X 10"'^ 

0 

2720.1 

48.78 

8000 

0.15 

D 

9.93 X 10"'^ 

15 

2720.1 

48.78 

8000 

0 

E 

9.93 X 10"'^ 

5 

2720.1 

48.78 

8000 

0 

F 

3.56 X 10"'^ 

0 

1986.3 

48.10 

8000 

0.05 

G 

3.56 X 10"'^ 

10 

1986.3 

48.10 

8000 

0 


Table 1 The model parameters in model A-G. 


3.1 Model A 

In Model A, we simulate a stellar bow shock in a uniform medium with the number density of uq = 
8000 cm~^. The velocity of the moving star is assumed to be u* = 10 km s~^. When simulating bow 
shock models, we carry out the calculation in the rest frame of the star first, so that the same procedures can 
also be used in the champagne flow cases. After the calculation, we convert all velocities to the values in 
the rest frame of molecular clouds. The velocities presented in this paper are all in the frame of reference 
of molecular clouds. In all models, the z-axis is parallel to the symmetrical axis and the positive direction 
is from the tail to the head of the cometary region. The star is at the position of {x,z) = (0,0). The 
effective temperature of the star is 40, OOOAT. The numbers of ionization photons (hv > 13.6eU) and the 
photodissociation photons (11.26 < hv < 13.6eU) emitted from the star per second are and 

1048 . 76 g-i, respectively. The mass-loss rate is M = 9.93 x lO“^M 0 ?/r“^, and the terminal velocity of the 
stellar wind is Vw = 2720.1 km These parameters are consistent with a star of mass M* = 40.9 Mq 
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Our simulation is stopped at 120,000 yr. Before that time, the ionization front ahead of the star has 
been approximately motionless relative to the star for a few 10^ years. In Figure 1, the number density of 
all materials, iJ+, and ions in model A are presented. A stellar wind bubble of low density 

(n < 5 cm~^) can be seen around the star. The bubble is surrounded by the photoionizated region of high 
density (n ~ 200 — 20000 cm~^), and a dense neutral shell (n ~ 10^“® cm~^). In the right panel of 
Figure 1, the density distributions of A^e"*" and ions are shown in grey scales, the Ne~^ and 

ions distribute in the whole H 11 region, but we only show the number density higher than 10“^ ® cm“® 
in order to highlight the different distributions between Ne~^ and ions. Although the high densities 

in the head of the H II region, which is defined as the ionized region of 2 > 0, leads to the high densities 
of both A^e’*' and there, the density of A^e"*" ions roughly increases with the distance from the star 

in the head of the H II region. The density of by contrast, decreases with the distance in the head 

region. Large fluxes of the photons with energy higher than the ionization potential of Ne~^ can increase 
the ionization rate from Ne~^ into Ne^'^ at the locations near the star. On the contrary, high electron density 
in the head region would increase the recombination rate. This leads to the different density distributions of 
A^e"*" and Ne^'^ ions. 



-1.5 -1 -0.5 0 0.5 -1.5 -1 -0.5 0 0.5 


Z [pc] Z [pc] 

Fig.l Gas density at logarithmic scales at the age of 120, 000 yrs. In the left panel, the top half 
of the figure shows the total density of the gas and the bottom half shows the density of ions. 
The right panel shows the densities of Ne^ ions (top half) and TVe^"*" ions (bottom half). 


In Figure 2, the velocity field of photoionized gas is shown. The result of the bow shock (model A) is 
shown in the left panel. A large proportion of ionized gas in the head of the H II region has red-shifted 
velocities at the direction from the tail to the head, and the velocities of the ionized gas in the boundary 
region, which is defined as a 0.025pc wide layer in the H II region near the ionization front, are approx¬ 
imately perpendicular to the ionization front. Particularly, the gas in front of the star moves at the similar 
direction as the stellar motion. The directions of motions of the ionized gas in the rest part of the HII region 
are gradually turned toward the tail. This phenomenon is because of the accelerations toward the tail due to 
the pressure density in the FI II region. The details of this reason is discussed in Section3.3.2. 

In Figure 3, the profiles of the [Ne II] line, the H30a line and the [Ne III] line from the H II region for 
different inclinations (9 = 0°, 30°, and 60° from the z-axis ) are presented. The profiles are all asymmet- 
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Z [pc] Z [pc] 

Fig. 2 The velocity fields of photoionized gas in model A (bow shock, left) and model B (cham¬ 
pagne flow, right). The velocity fields in the stellar wind bubble(n < 50 cm~^) and outside the 
H II region (X(iJ+) < 0.1) are not shown. The arrows represent velocities of u > 1 km s~^, 
and the lengths of them are proportional to their absolute values. The position of the star is (0,0). 
The five contour levels are at 50, 500, 2000, 20000 and 100000 cm~^. 



A V [km/s] 



A V [km/s] 



A V [km/s] 


Fig. 3 profiles of the [Ne II] 12.81/rm line (left), the H30a line (middle) and the [Ne III] 
15.55p,m line (right) from the cometary H II regions in model A for three inclination angles 
(0°, 30°, 60°). 


Line 

Inclination 

Peak 

(km s~^) 

FWHM 

{km s~^) 

FWCV 

{km s~^) 

Luminosity 
{erg s"!) 

[Ne II] 12.81^J,m 

0° 

3.1 

20.8 

-3.32 

7.75 X 10^® 


CO 

o 

0 

2.5 

19.0 

-2.92 



60° 

0.9 

13.0 

-1.70 


H30a 

0° 

-1.4 

36.8 

-3.40 

1.07 X 10^° 


CO 

o 

0 

-1.4 

35.6 

-2.99 



o 

0 

-I 

33.2 

-1.74 


[Ne III] 15.55/im 

0° 

1.7 

17.6 

-3.28 

1.03 X 10^® 


CO 

o 

0 

1.7 

16.6 

-2.85 



60° 

0.9 

13.0 

-1.64 



Table 2 Peaks, FWHMs and flux weighted central velocities(FWCV) of lines at different angles 
in model A 
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Fig. 4 the line profiles of the [Ne II] 12.81^m from the whole HII region, the head and boundary 
region and the rest part of the H II region. The left panel is for model A, and the right panel is 
for model B. 


profiles are skewed to the right and have a long tail at the left side. With increasing inclination angles, these 
line profiles become less asymmetrical and narrower. It is obvious that the H30a line is much broader than 
the [Ne II] line and the [Ne III] line. This is due to the larger thermal broadening of lower mass hydrogen. 

In Table 2, the peak locations, FWHMs and flux weighted central velocities (FWCV) of the lines for 
three inclinations of 0°, 30° and 60° are shown. For the [Ne II] line, the peak locations and FWHMs of 
the profiles decrease with increasing inclination angles. Although the peak locations are all red-shifted 
(3.1, 2.5 and 0.9 km s~^), the flux weighted central velocities are all blue-shifted (—3.32, 2.92 and — 
1.70 km s~^). This suggests that the contributions of blue-shifted tail in the profile are significant. Because 
the line profile for the inclination of 0° seem to be consist of two components, we calculate the [Ne II] 
line profile from the head and the boundary region and the line profile from the rest part of the H II region, 
separately. These two line profiles are plotted in Figure 4. It can be seen that the emission from the head 
region and the boundary region leads to the red-shifted peak location of the total line profile, and and the 
gas from the rest part of the H II region mainly contributes to the blue-shifted tail of the line profile. The 
luminosity of the [Ne II] line is 7.75 x 10^^ erg s~^ which is independent of inclinations. 

For the [Ne III] line, the peak locations are slightly red-shifted (1.7, 1.7 and 0.9 km s“^), and the flux 
weighted central velocities are all blue-shifted (—3.4, — 2.99 and — 1.74 km s“^). The reason of this is 
the same as in the [Ne II] line profile. The line luminosity of the [Ne III] line is 1.03 x 10^® erg This 
value is larger than that of the [Ne II] line. For the H30a line, both the peak locations (—1.4, — 1.4 and — 
1.0 km and the flux weighted central velocities (—3.28, — 2.85 and — 1.64 km s“^) are blue-shifted. 
This is due to the large thermal broadening of lower mass hydrogen. By increasing the broadening, the peak 
locations of all lines will approach the flux weighted central velocities. The blue-shift of the flux weighted 
central velocities of the hydrogen line suggests the blue-shifted gas motion dominates for the ionized gas 
in bow shock model. It seems odd that the ionized gas moves mainly at a direction opposite to the direction 
of the star and the shock structure. The reason is also the accelerations toward the tail due to the pressure 
density in the H II region. The luminosity of the H30a line is 1.07 x 10^° erg s~^ in model A. This is 
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3.2 Model B 

In model B, the evolution of a champagne flow including a stellar wind is simulated. The density distribution 
follows an exponential law as n{z) = noexp{z/H). z is the axial coordinate along the symmetrical axis 
and towards the center of molecular cloud. The density at the position of the motionless massive star (z = 0) 
is no = 8000 cm~^ initially and the scale height is H = 0.05 pc. The parameters of the massive star and 
the stellar wind are the same as in model A. In model B, the simulation of the time evolution is ceased at 
160,000 yr when the line profiles are roughly stable. The H II region and the neutral region should reach 
approximate pressure equilibrium, and the champagne flow has completely cleared the low-density material 
from the grid at the time. 

The density distribution in model B at the end of simulation is presented in Figure 5. As in model 
A, there are also a stellar bubble surrounded by the photoionized region and a dense neutral shell which 
separates the H II region from the dense cloud. The size of H II region in model B is bigger than in model 
A, but the total number of ions is only 0.74 of the number in model A. Hence, the average density of 
electrons in model B is lower than in model A. The low electron density decreases the recombination rate 
of Ne'^~^ so that the number of the Ne~^ ions is only 0.38 of the number of the ions in model B. In 
the tail, for the same reason, the relative fraction of Ne^ is generally lower than 0.3 in model B while that 
is always higher than 0.6 there in model A. The difference in the numbers between the Ne^ and Ne^'^ 
ions is easily found in the right panel of Figure 5 where the density distributions of and ions 

are presented. The velocity field of photoionized gas in model B is presented in the right panel of Figure 2. 
In the boundary region, the axial components of the velocities are low, and the ionized materials in the rest 
part of the H II region all have an apparently blue-shifted velocity for the inclination of 0°. 

The line profiles from the H II region are plotted in Figure 6. In model B, the line profiles are also 
asymmetrical as in model A. For the [Ne II] line, we separate the boundary region from the rest part of the 
H II region because the axial velocity gradually change along the symmetrical axis from ^ 0 km at 
the boundary to ~ —20 km at the tail. The line profiles from these two regions are plotted in the right 
panel of Figure 4. The [Ne II] line profile is consist of two components; a narrow component with a slightly 
blue-shifted peak location contributed by the ionized gas in the boundary region and a highly blue-shifted 
broad component with the peak location at z; = —9.7 km contributed by gas from the rest part. This 
suggests that gas in boundary region is not accelerated too much. Meanwhile, the ionized gas in the rest part 
has been accelerated to high velocities toward the tail direction. When the inclination angle increases, the 
highly blue-shifted component is affected strongly and the center of the component moves toward the red 
but a little change happens to the other component. It is obvious that the profiles of the [Ne II] line in model 
B are more biased to the left side than those in model A. In addition, because the number density is smaller 
than the critical densities of [Ne II] and [Ne III] lines in the H II region, the collisional de-excitation is not 
important for these lines. The fluxes of lines from a unit volume is fui{i) oc rig. So the number of the [Ne II] 
line photons emitted from the dense boundary region is not negligible relative to the total number although 
the number of the Ne^ in the boundary region is just 0.08 of the total number. For the H30a line and the 
[Ne III] line, the proportions of 77+ and A^e^+ ions in the boundary region are smaller. So the proportions 
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In Table 3, the peak locations, the FWHMs and the FWCVs in model B are provided. For the [Ne II] line 
profiles, the peak locations (—8.9, —6.3, —2.5 fcm s“^) and the FWCVs (—8.52, —7.38, —4.28fcms“^) 
are all blue-shifted. The values of the peak locations and the FWCVs in model B become less blue-shifted 
with increasing inclination angles, and they are all much more blue-shifted than the corresponding values 
in model A. For the H3Qa recombination line, we find that both the peak locations (—9.4, — 7.8, — 
4.2 km and the FWCVs (—10.13, — 8.77, — 5.06 km s~^) are much more blue-shifted than those 
in model A because the main part of the ionized gas in model B are forced to flow away from the molecular 
cloud due to the density gradient. For the [Ne III] line, both of the FWCVs (—10.81, —9.37, —5.45 km s“^) 
and the peak locations (—9.3, — 7.5, — 3.7 km s~^) of the profiles are easily distinguished from those 
in model A. This suggests the ionized gas in the inner part of the H II region has a high velocity (~ 
— 10 km s~^) towards the tail. Also in Table 3, the [Ne III] line profiles have the most blue-shifted peak 
locations and FWCVs in the three lines. We find that the three line luminosities are lower in model B than 
in model A. This is due to the lower density in model B. 



Fig. 5 Gas density at logarithmic scales at the age of 160, 000 yr in model B. The top half in the 
left panel shows the density of all materials. The bottom half in the left panel shows the density 
of 77+ ions. The top half in the right panel shows the density of A^e+ ions. The bottom half in 
the right panel is for 7Ve^+ ions. The density is in the units of cm~^. 



A V [km/s] 



A V [km/s] 



A V [km/s] 


Fig. 6 profiles of the [Ne II] 12.81^m line (left), the 7730a line (middle) and the [Ne III] 
15.55/iTO line (right) from the cometary H II regions in model B for three inclination angles 
(0°, 30°, 60°). 
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3.3 Comparison between Champagne Flow and Bow Shock Models 

We have calculated model A as a bow shock model and model B as a champagne flow model and obtain 
the line profiles, peak locations and FWCVs of the [Ne II] 12.81^m line, the H30a line and the [Ne III] 
15.55pm line. It is worth to point out that line profiles and peak locations can be influenced by broadening. 
For example, if we change the broadening to be 5 km s~^ , the peak location of the [Ne II] line profile for 
the inclination of 0° in model A will be changed to be 1.1 km s~^ from 3.1 km But the FWCVs are 
independent of the broadening. 

For the H30a line, both of the peak locations and the FWCVs in model A are much less blue-shifted 
than in model B for every inclinations. The FWCVs in model A are all higher than —3.3 km s~^, but those 
in model B are all lower than —5 km s~^. For the [Ne III] line, the difference is more apparent. The peak 
locations in model A are red-shifted while those in model B are blue-shifted for every inclinations. The 
FWCVs of the [Ne III] line in model A are all higher than —3.5 km s~^, but in model B they are all lower 
than —5.4 km s~^. Especially for the inclination of zero, the difference of the FWCVs in these two models 
is the largest and equal to 7.41 km s~^. For the [Ne II] line, it is also the case that the emission is generally 
more blue-shifted in model B than in model A. And the FWCVs and the line profiles between these two 
models are also much different. 

So it is possible to distinguish a champagne flow from a bow shock by using the [Ne II] line. And the 
difference between the champagne flow and bow shock is more obvious in the H30a line and the [Ne III] 
line. In order to check whether these conclusions are applicable in more general cases of champagne flows 
and bow shocks, we compute other models and test the effects of the density gradient, the velocity of the 
moving star, the mass of the star on the line profiles of the three lines. 

3.3.1 Density Gradient in the Champagne Flow 

A Champagne flow model with a shallow density gradient (H = 0.15pc) is computed in model C. The other 
parameters are kept same as in model B. We cease the evolution in model C at 160, OOOyr as in model B. 
The line profiles are also roughly stable, and the pressure equilibrium has been formed. 

The line profiles and properties of the H II regions for model B and C are shown in Figure 7 and Table 
3. In model C, the H II regions can also be divided into a boundary region and the rest region as in model 
B. The velocities in the boundary region are only slightly blue-shifted and close to zero, and the gas in the 
rest part of the H II region has a obviously blue-shifted velocity. Because of the shallow density gradient, 
the peak location of the [Ne II] line profile from the rest part of the H II region is at —7.2 km s“^. This 
suggests that acceleration of the ionized gas is smaller in model C than in model B. Beside this, the line 
profiles for model C are similar to those for model B. 

In Table 3, we find that the peak locations and the FWCVs of the three lines are all blue-shifted. The 
values of the FWCVs and the peak locations in model C are slightly less blue-shifted than the corresponding 
values in model B. And as in model B, the FWCVs and the peak locations of the [Ne III] line have the most 
blue-shifted values in the line profiles of the three lines. 

Our calculations suggest that the line profiles become more blue-shifted with the increasing density 
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Fig. 7 The profiles of the [Ne II] 12.81fim line (left panels), the HiOa line (middle panels) 
and the [Ne III] 15.55pm line (right panels) from the H II regions for three inclination angles in 
model B (top panels) and model C (bottom panels). 


Model BH = O.Obpc 


Line 

Inclination 

Peak 

FWHM 

FWCV 

Luminosity 



{km s~^) 

{km s~^) 

{km s~^) 

(erg s”!) 

[Ne II] 12.81^™ 

0° 

-8.9 

15.4 

-8.52 

3.43 X 10^5 


CO 

o 

0 

-6.3 

13.8 

-7.38 



60° 

-2.5 

10.8 

-4.28 


H30a 

0° 

-9.4 

36.0 

-10.13 

6.08 X 10^® 


CO 

o 

0 

-7.8 

35.6 

-8.77 



60° 

-4.2 

35.6 

-5.06 


[Ne III] 15.55/im 

0° 

-9.3 

11.4 

-10.81 

6.60 X 10®® 


CO 

o 

0 

-7.5 

11.8 

-9.37 



60° 

-3.7 

12.0 

-5.45 


Model EH = 0.15pc 


Line 

Inclination 

Peak 

FWHM 

FWCV 

Luminosity 



{km s~^) 

{km s~^) 

{km s~^) 

{erg 

[Ne II] 12.81tJ,m 

0° 

-6.3 

15.6 

-6.60 

3.68 X 10®® 


CO 

o 

0 

-4.7 

13.6 

-5.70 



0 

o 

CD 

-2.1 

9.8 

-3.29 


H30a 

0° 

-1.4 

36.4 

-8.32 

4.90 X 10®® 


CO 

o 

0 

-6.2 

36.0 

-7.21 



60° 

-3.8 

34.4 

-4.16 


[Ne III] 15.55/im 

0° 

-6.7 

12.2 

-9.48 

4.70 X 10®® 


CO 

o 

0 

-5.9 

11.2 

-8.24 



60° 

-2.9 

10.2 

-4.77 



Table 3 Peaks, FWHMs and flux weighted central velocities(FWCV) of lines at different angles 
in model B and C. 
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3.3.2 Stellar Motion in the Bow Shock 


The difference among model A, D, and E is the stellar velocity with respect to the ambient molecular clouds 
(u* = 10, 15 and 5 km s~^). The line profiles are plotted in Figure 8. The peak locations and the FWHMs 
are shown in the Table 4. We find that the line profiles are more biased towards the red-shifted side with the 
increasing stellar velocity. The peak locations and the FWCVs of the three line profiles in model D are all 
red-shifted. On the contrary, these values in model E are all blue-shifted. This is because the proportion of 
ionized gas with a red-shifted velocity is the highest in model D and is the lowest in model E. As in model 
A, the red-shifted materials are mainly in the boundary region and the head region. In model D, the [Ne II] 
line profile is also consist of the components from the red-shifted region and from the blue-shifted region. 
But this is not obvious in the profiles of the [Ne III] and the H30a line due to the different relative fraction 
of and and the large broadening of iL+. In model E, because the proportion and the velocities of 
red-shifted ionized gas are both low, the contribution to the red-shifted part of the [Ne II] line profile is small 
as well. In our results of the simulations, the size of the H II region decreases with the increasing stellar 
velocity. A shorter distance due to the smaller size causes the flux of the ionizing photons to be stronger so 
that the ionization rate from Ne~^ to rises. But the density of electrons could be higher and increase 

the recombination rate of Ne^~^ ions. The comparison of the line luminosities of the [Ne II] line and the 
[Ne III] line shows that higher stellar velocity causes lower ratio of the [Ne II] line luminosity to the [Ne 
III] line luminosity. This suggests that the increase of the ionization rate from Ne~^ to is more than 

that of the recombination rate. 


In bow shock models, the ionized gas compresses the neutral and cold materials ahead of the moving 
star into a dense shell while the density of the ionized gas in HII region is much lower than the density of the 
shell. When the evolution reaches quasi-steady state, the ionization front and the shock front are motionless 
relative to the star. Along the arched shell from the apex to the tail, the axial velocity of the neutral gas 
decreases from the stellar velocity to zero. In the ionized region, the gases just peeled from the shell by the 
ionization are dense 10000 cm~^) and have a velocity slightly less than the velocity of nearby neutral 
gases in the shell. With the expansion and moving to the tail of the H II region, the densities of the ionized 
gases gradually decrease to lower values 50 cm~^). This leads to a pressure gradient from the head to 
the tail as in a champagne flow model. Although the ionized gases ahead of the star are accelerated by the 
stellar wind, the advection in the head region pushes these gases to the sides where the pressure gradient 
is dominant. So if the stellar velocity are higher, it need more time for the pressure gradient to accelerate 
the ionized gas from shell toward the blue-shifted direction, so that the proportion of the red-shifted ionized 
gas in the H II region will be higher. In model E, the effect of the pressure gradient is dominant in the H II 
region due to the low stellar velocity. Hence, the line profiles for model E in Figure 8 are more biased to the 
left side than those in model A and model D. And the peak locations and the FWCVs in model E are not 
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Fig. 8 Same as Figure3, but for model D(top panels) and model E(bottom panels). There are the 
profiles of the [Ne II] 12.81/im line (left panels), the H30a line (middle panels) and the [Ne 
III] 15.55pm line (right panels) from the cometary H II regions in model A for three inclination 
angles (0°, 30°, 60°). 


3.3.3 the Mass of the star 

We assume a less massive star (M* = 21.9 Mq) in model F and G. This leads to a low effective temperature 
35000AT, a weak stellar wind and a weak radiation (See Table 1 ). The other parameters in model F and model 
G are kept same as in model B and model A, respectively. 

In model F, because of the weaker stellar wind and ionizing radiation, the size of the H II region is 
also smaller than in model B. The lower effective temperature causes the relative fraction of ions to 
be higher. For example, in the tail, X{Ne^) is generally higher than 0.6 in model F but lower than 0.3 in 
model B. 

The H II region can also be divided into the boundary region and the rest region. For the [Ne II] line, 
the proportion of the photons emitted from the boundary region is lower than in model B. This is due to the 
higher relative fraction of A^e+ in the inner part of the H II region and causes the lower normalized flux at 
the velocities close to 0 in the [Ne II] line profile. The other two line profiles in model F are similar to those 
in model B. 

In Table 5, the properties of the line profiles are shown. The FWCVs in model F are more blue-shifted 
than the corresponding values in model B. This is related with the higher relative fraction of Ne^ in the 
inner part and the lower fluxes emitted from the boundary region where the velocities close to zero. In 
addition, the low effective temperature and the ionization radiation of the star lead to similar distributions 
of Ne^ and ions in the H II region so that the FWCVs of the [Ne II] line, the H30a line and the [Ne 
III] line are closer to each others in model F than in model B. Because of the high relative fraction of Ne~^, 
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Model D = 15/cm s ^ 


Line 

Inclination 

Peak 

FWHM 

FWCV 

Luminosity 



{km s~^) 

{km s~^) 

{km s~^) 

(erg 

[Ne II] 12.81/.tm 

0° 

4.1 

25.2 

0.73 

7.79 X 10^® 


30° 

3.5 

19.2 

0.63 



60° 

2.7 

13.6 

0.36 


H30a 

0° 

5.8 

36.4 

1.14 

1.18 X 10®° 


30° 

3.8 

35.6 

0.98 



60° 

1.0 

34.4 

0.57 


[Ne III] 15.55^^m 

0° 

3.9 

18.8 

1.26 

1.14 X 10®° 


30° 

3.5 

16.4 

1.08 



60° 

2.7 

13.4 

0.62 


Model E u* = 5fcm s~^ 


Line 

Inclination 

Peak 

FWHM 

FWCV 

Luminosity 



(km s~^) 

{km s~^) 

{km s~^) 

(erg s”®) 

[Ne II] 12.8lAtm 

0° 

-1.9 

13.2 

-5.49 

8.24 X 10®® 


30° 

-1.9 

12.6 

-4.80 



0 

o 

CC 

-1.3 

10.6 

-2.80 


H30a 

0° 

-3.8 

36.4 

-5.09 

9.54 X 10^° 


CO 

o 

0 

-3.4 

36.0 

-4.41 



60° 

-2.2 

34.4 

-2.55 


[Ne III] 15.55^^m 

0° 

-2.3 

11.4 

-4.78 

8.09 X 10®® 


30° 

-1.9 

11.2 

-4.17 



60° 

-1.3 

10.4 

-2.43 



Table 4 Peaks, FWHMs and flux weighted central velocities (FWCV) of lines in model D and 
E. 


The weaker EUV flux and the weaker ram-pressure of the stellar wind in model G result in a smaller H 
II region and result in a smaller mass ratio between the head region and the entire H II region than in model 
A (0.11/0.15 in model G/A). The proportion of the red-shifted ionized gas is also lower than in model A 
(0.13/0.18 in model G/A). The lower effective temperature leads to a higher relative fraction of Ne^ ions 
in the H II region as in model E. 

The peak locations, line luminosities and the flux weighted central velocities in model G are also pre¬ 
sented in Table 5. The values of the EWCVs of the three lines are close to each others for every inclination 
as in model E. But it is only a coincidence that the EWCVs of the [Ne II] line are approximately equal to 
those of the [Ne III] line. The EWCVs in model G are all more blue-shifted than the corresponding values 
in model A. This is attributed to the lower proportion of the red-shifted ionized gas mentioned above. It is 
also obvious that the line luminosity of the [Ne II] line presented is much stronger than that of the [Ne III] 
line in model G due to the low effective temperature. 

3.4 Line Profiles Computed by Using a Slit 

We have calculated four bow shock models, but only the EWCVs and the peak locations in model D with a 
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Fig.9 There are the profiles of the [Ne II] 12.81^m line (left panels), the HiOa line (middle 
panels) and the [Ne III] 15.55/iTO line (right panels) from the cometary H II regions for model F 
(top panels) and model G (bottom panels). 


model F M* = 21.9Mq H 

= O.OSpc 


Line 


Inclination 

Peak 

FWHM 

FWCV 

Luminosity 




{km s~^) 

{km s~^) 

{km s~^) 

(erg s”!) 

[Ne II] 12.81/im 


0° 

-9.3 

17.0 

-10.48 

6.07 X lO^'^ 



CO 

o 

0 

-7.7 

15.2 

-9.10 




60° 

-2.5 

II.6 

-5.27 


H30a 


0° 

-10.2 

37.2 

-II.OI 

5.59 X 10^® 



CO 

o 

0 

-8.6 

36.4 

-9.54 




60° 

-5.0 

35.2 

-5.51 


[Ne III] 15.55/im 


0° 

-8.5 

12.6 

-11.33 

3.36 X 10®^ 



CO 

o 

0 

-7.1 

12.4 

-9.86 




60° 

-3.1 

II.2 

-5.73 


model G M* = 21.9M0 w* 

= 10km s~^ 


Line 


Inclination 

Peak 

FWHM 

FWCV 

Luminosity 




{km s~^) 

{km s~^) 

{km s~^) 

{erg 5“^) 

[Ne II] 12.81/im 


0° 

4.3 

21.0 

-3.82 

2.92 X 10®® 



CO 

o 

0 

2.3 

18.4 

-3.34 




0 

o 

CD 

0.5 

12.8 

-1.93 


H30a 


0° 

-1.4 

38.8 

-3.59 

2.25 X 10^9 



CO 

o 

0 

-1.4 

37.6 

-3.II 




60° 

-I.O 

35.2 

-1.80 


[Ne III] 15.55/im 


0° 

O.I 

18.4 

-3.82 

9.46 X lO®"' 



CO 

o 

0 

0.7 

16.8 

-3.34 




60° 

0.3 

13.6 

-1.94 



Table 5 Details of the line profiles in model F and model G. 
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because the velocities in the low-density and large-volume inner part of the H II region are always blue- 
shifted. As is mentioned in Sections. 3.2, the ionized gas just ionized from the apex of the ached shell has 
a similar velocity to the stellar velocity and is dense. So if we compute the line profiles from a slit along 
the symmetrical axis of the projected 2D image rather than from the whole H II region, the influence of 
the dense gases in the head of the H II region will be highlighted while the emission from the low-density 
region will be weakened. The method is applied for all the models, and the profiles of the [Ne II] line are 
presented in Figure 10. The peak locations of the profiles of the [Ne II] line and the H30a line are shown 
in Table 6. 

For the [Ne II] line, the profiles are all single-peaked. In bow shock models (model A, D, E and G), 
the peak locations are slightly less than cos(0)u*. In champagne flow models (model B, C and F), it is 
different that the peak locations of the [Ne II] line profiles are close to zero. The red-shifted peak locations 
in model C suggest that the propagation velocity of the ionization front in model C is a little faster than in 
other champagne flow models at the age of 160,000 yr. For the H30a line, because of the large thermal 
broadening, the relation of the peak locations to the stellar velocities in bow shock models is less explicit. 
But, the champagne flow models are easily distinguished from the bow shock models due to the red-shifted 
peak locations in bow shock models and the blue-shifted values in champagne flow models. 

When comparing the FWCVs presented in Table 6, we find that the FWCVs in champagne flow models 
are all blue-shifted. In bow shock models, the FWCVs are obviously red-shifted when the stellar velocity is 
high. But, the FWCVs in model E are slightly blue-shifted due to the low stellar velocity (v^, = 5 km s~^). 
And it suggests that the EWCVs in a bow shock model with a stellar velocity lower than 5 km s~^ are not 
easily distinguished from the corresponding values in champagne flow models. 



Fig. 10 There are the profiles of the [Ne II] 12.81p,m line from the slit along the symmetrical 
axis of the projected 2D image of the models for three inclinations. The top panels are plotted 
for the bow shock models, and the bottom panels are for the champagne flow models. 


3.5 the Flux Weighted Central Velocities 


As is mentioned in Sections. 2, the number density is n ~ rie <C Ucr- And in our calculation, the temperature 

























































Peak(/cm s“^) 

[Ne 11] 

12.81/Tm 


H30a 

Inclination 

0° 

30° 

60° 

0° 

CO 

O 

0 

60° 

Model A 

8.6 

7.0 

3.1 

7.8 

6.6 

3.8 

Model B 

-0.1 

-0.2 

-0.5 

-2.6 

-2.2 

-1.4 

Model C 

1.4 

0.3 

0.3 

-1.4 

-1.4 

-0.6 

Model D 

12.6 

12.0 

7.0 

8.6 

7.4 

4.2 

Model E 

4.8 

4.1 

2.0 

2.6 

2.2 

1.0 

Model F 

-0.9 

-1.2 

-1.1 

-3.0 

-2.6 

-1.4 

Model G 

8.5 

7.4 

4.2 

lA 

6.2 

3.4 

FWCV(A;ms-l) 

[Ne 11] 

12.81//m 


H30a 

Inclination 

0° 

30° 

60° 

0° 

CO 

O 

0 

60° 

Model A 

4.18 

3.50 

2.07 

4.92 

4.26 

2.46 

Model B 

-3.82 

-3.24 

-1.77 

-4.95 

-4.29 

-2.47 

Model C 

-2.42 

-2.09 

-1.24 

-4.29 

-3.72 

-2.15 

Model D 

7.27 

6.46 

3.86 

6.14 

5.32 

3.07 

Model E 

-0.12 

-0.30 

-0.57 

-1.00 

-0.86 

-0.50 

Model F 

-5.38 

-4.60 

-2.60 

-5.70 

-4.94 

-2.85 

Model G 

6.08 

5.18 

2.89 

5.16 

4.47 

2.58 


Table 6 The peak locations and the flux weighted central velocities of the [Ne II] 12.81/rm line 
and the H30a line from the slit along the symmetrical axis of the projected 2D image. 


photoionized region, the flux of the [Ne II] and the [Ne III] line from a unit volume is approximately as 
follow: 


bkNkAuihvui for H30a, k = 30, i = 0 

fui(i) = { AbNeX(Ne^){nl/{2ncr))Auihiyui /or [Ne II], i = 1 ■ (4) 

AbNeX{Ne^~''){5nl/{3ncri))Ay,ihvui for [Ne III], i = 2 
Since the neon atoms are almost all ionized in the H II region, we can get the relation that X{Ne~^) + 
X{Ne‘^^) = 1. Other parameters in Eq. (4) are constant due to the constant temperature in photoionized 
region. Hence, the fluxes from a unit volume of the three lines are all proportional to rig. Then we can obtain 
the following relation: 


(: 


L(l) 


r^(i) 


L(2) 


r^(2))/(. 


L(l) 


L(2) 


-) = ^(0) , (5) 


'1.302 x 10-21 1.162x10-21 '1.302 X 10-21 1.162 x 10-21' 

where L{1) and L{2) are the line luminosities of the [Ne II] line and the [Ne III] line, respectively. 
u(0), i;(l) and i;(2) are the FWCVs of the profiles of the H30a line, the [Ne II] line and the [Ne III] line, 
respectively. If the abundance of neon has been obtained, we can approximately compute the line luminosity 
and the FWCVs of the [Ne III] line from those of the other two lines: 


F(l) 


Ne 


L{2) = (L(0) 


1.801 X lO^^Ab 


) X 1.608 X 10^°AbNe , 


( 6 ) 
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where L(0) is the line luminosity of the H30a line. The [Ne III] line is not detectable from the ground, 
but the [Ne II] line and the H30a line can be observed with the ground based telescopes. These equations 
could be useful since the difference in the [Ne III] line between the bow shock models and the champagne 
flow models is the most obvious among the three lines. For observing the [Ne III] line, a typical resolution 
(A/AA = 30, 000) is needed to resolve the difference in the [Ne III] line profiles. The line also needs to be 
observed by using space telescopes since it can not be observed from the ground. In our models, the average 
deviation of the approximate values computed with the Eq. (6) and (7) from the accurate values is 6%. 

4 CONCLUSIONS 

When considering the gas kinematics of compact H II region models, people generally predict that the line 
emission from the ionized region may be red-shifted along the line of sight for a bow shock and blue-shifted 
for a champagne flow when looking into the regions from the tail. Our simulation proves that this envision 
is generally true, but it is violated in the case of low/medium stellar speed bow shock. When the stellar 
velocity is lower than 5 km it is difficult to distinguish the bow shock from a champagne flow. On the 
contrary, A bow shock is very easily distinguished from champagne flows if the stellar velocity is higher 
than 10 km s~^. 

In this paper, we have simulated the evolution of the cometary H II region by bow shock models and 
champagne flow models, and have displayed the density distribution and velocity fields in bow shock models 
and champagne flow models. We have studied the [Ne II] 12.81/rm line, the H30a line and the [Ne III] 
15.55^m line profiles from these models. We make comparisons of the line profiles in bow shock models 
and champagne flow models. We And that: 

1. After the champagne models reach quasi-steady state, the line profiles are all biased towards the 
blue-shifted side and have obviously blue-shifted peak locations and flux weighted central velocities. The 
peak locations and the flux weighted central velocities of the [Ne II] line profiles are the least blue-shifted 
in the three lines meanwhile those of the [Ne III] line profiles are the most blue-shifted in champagne flow 
models. In our champagne flow models, when the inclination angle is no more than 60°, the flux weighted 
central velocities are lower than —3.0 km s~^ for the [Ne II] line, are lower than —4.0 km s~^ for the 
H30a line and are lower than —4.5 km s~^ for the [Ne III] line. 

2. For bow shock models. We And that the line profiles become more biased towards the red-shifted side 
with increasing stellar velocity. Except for the model D with a star of u* = 5 km s~^, in the bow shock 
models computed in this paper, the flux weighted central velocities are higher than —3.82 km s~^ for the 
[Ne II] line and the [Ne III] line, and are higher than —3.6 km s~^ for the H30a line. 

3. The peak locations of the [Ne II] line profiles from the slits along the symmetrical axis of the projected 
2D image are slightly less than cos(0)u* in bow shock models and no more than 1.5 km s~^ in champagne 
flow models. The peak locations of the H30a line profiles from the slits are less red-shifted or more blue- 
shifted than those of the [Ne II] line. 

We have shown the profiles of the [Ne II] line at 12.81fim, the H30a recombination line and the [Ne 
III] line at 15.55^m in this paper. These results are useful to recognizing the causes of the cometary H II 
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